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Abstract 

Background: Chitosan oligosaccharide (COS), a deacetylated derivative of chitin, is an abundant, and renewable 
natural polymer. COS has higher antimicrobial properties than chitosan and is presumed to act by disrupting/ 
permeabilizing the cell membranes of bacteria, yeast and fungi. COS is relatively non-toxic to mammals. By 
identifying the molecular and genetic targets of COS, we hope to gain a better understanding of the antifungal 
mode of action of COS. 

Results: Three different chemogenomic fitness assays, haploinsufficiency (HIP), homozygous deletion (HOP), and 
multicopy suppression (MSP) profiling were combined with a transcriptomic analysis to gain insight in to the mode 
of action and mechanisms of resistance to chitosan oligosaccharides. The fitness assays identified 39 yeast deletion 
strains sensitive to COS and 21 suppressors of COS sensitivity. The genes identified are involved in processes such 
as RNA biology (transcription, translation and regulatory mechanisms), membrane functions (e.g. signalling, 
transport and targeting), membrane structural components, cell division, and proteasome processes. The 
transcriptomes of control wild type and 5 suppressor strains overexpressing ARL1, BCK2, ERG24, MSG5, or RBA50, 
were analyzed in the presence and absence of COS. Some of the up-regulated transcripts in the suppressor 
overexpressing strains exposed to COS included genes involved in transcription, cell cycle, stress response and the 
Ras signal transduction pathway. Down-regulated transcripts included those encoding protein folding components 
and respiratory chain proteins. The COS-induced transcriptional response is distinct from previously described 
environmental stress responses (i.e. thermal, salt, osmotic and oxidative stress) and pre-treatment with these well 
characterized environmental stressors provided little or any resistance to COS. 

Conclusions: Overexpression of the ARL1 gene, a member of the Ras superfamily that regulates membrane 
trafficking, provides protection against COS-induced cell membrane permeability and damage. We found that the 
ARL1 COS-resistant over-expression strain was as sensitive to Amphotericin B, Fluconazole and Terbinafine as the 
wild type cells and that when COS and Fluconazole are used in combination they act in a synergistic fashion. 
The gene targets of COS identified in this study indicate that COS's mechanism of action is different from other 
commonly studied fungicides that target membranes, suggesting that COS may be an effective fungicide for 
drug-resistant fungal pathogens. 
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Background 

Chitin is an abundant natural polymer, second only in 
biomass to cellulose [1]. It is a common constituent 
of crustacean exoskeletons and arthropod cuticles [2] as 
well as the cell walls of most fungi [2,3] . There is an esti- 
mated 10 gigatons of chitin recycled in nature each 
year [4]. 

Chitosan is a polymer of N-glucosamine obtained by 
partial chitin N-deacetylation [1]. Acid hydrolysis or en- 
zymatic cleavage of the glycosidic linkages of chitosan 
chains [5], yields shorter chains of <1-10 KDa chitosan 
oligosaccharides (COS). These polymers are less than 
100 glucosamine monomers and are more water soluble 
and readily absorbed in vivo compared to chitosan [6] . 

Allan and Hadwiger [7] first showed the fungicidal ef- 
fect of chitosan, and since then, several studies have 
examined chitosan sensitivity in different fungi [8-13]. 
Chitosan inhibits both hyphal growth and spore germin- 
ation [14] and reduces toxin production by plant patho- 
genic fungi [15]. Recently, chitosan has been shown to 
increase conidiation in filamentous fungi at concentra- 
tions where hyphal growth was impaired [8]. Together 
these observations suggest that chitosan exerts its anti- 
fungal activity by multiple mechanisms. 

We used COS in our study rather than chitosan be- 
cause it is more soluble and biologically active. The bio- 
logical activity of chitosan and COS is dependent on its 
molecular weight, degree of deacetylation and pH of the 
medium [16]. High degrees of deacetylation and an 
acidic pH (when most amino groups are in the free 
base form) yield the highest activity against susceptible 
fungi [17]. 

The antibiotic activity of chitosan and COS [16,18-20] is 
likely due to the permeabilization of bacterial plasma 
membranes [18,21]. Plasma membrane damage has also 
been suggested as an explanation of the fungicidal effects 
of chitosan in yeast and filamentous fungi [10,12,22,23]. 
Plasma membrane permeabilization by chitosan has also 
been detected (indirectly) in Saccharomyces cerevisiae, 
where deletion of genes encoding for proteins involved in 
maintaining plasma membrane integrity increased chito- 
san sensitivity [23]. Although it has been speculated that 
plasma membrane permeabilization by chitosan in bac- 
teria and fungi is associated with the interaction between 
the positive amino groups of chitosan and the negative 
charges of phospholipids, no conclusive data are 
available [10,12,18]. 

To gain a better understanding of the mode of action 
of small (< 6 kDa) chitosan oligosaccharides (COS- 
5.44), and identify putative resistance mechanisms, we 
performed three chemogenomic screens [24,25]. These 
assays examined the effect of gene dosage to uncover 
bona fide targets of COS in yeast. Two chemogenomics 
assays interrogate the yeast deletion collection (~ 6000 



deletion strains) in a single culture in parallel. Each 
deletion strain in this collection contains a unique 20- 
base-pair DNA tag used to quantify the fitness of indi- 
vidual strains using an oligonucleotide array [24,26]. 
Both essential heterozygous and non-essential homozy- 
gous diploids were assayed. The haploinsufficiency pro- 
filing (HIP) assay identifies genes that show increased 
drug sensitivity in heterozygous deletion strains. The 
heterozygous deletion strain that is most sensitive to a 
given compound often indentifies its target [24,26]. The 
homozygous profiling (HOP) assay identifies genes that 
buffer the drug target pathway and are required for re- 
sistance to the compound [24,26]. A multicopy sup- 
pression profiling (MSP) screen, where genes are 
overexpressed to identify those that confer resistance to 
the compound of interest was also performed [25]. The 
genes identified from the chemogenomic assays were 
individually validated to confirm their involvement in 
COS sensitivity or resistance. For five of the strains 
that conferred resistance when overexpressed, we per- 
formed a transcriptional analysis. The ARL1 (ADP- 
ribosylation factor-like 1) gene was identified in the 
HIP-HOP and MSP chemogenomics assays as being 
involved in yeast resistance to COS-5.44. This gene 
encodes a soluble GTPase that is a member of the Ras 
superfamily and has been shown to be involved in 
membrane traffic regulation [27,28]. ARL1 overexpres- 
sion strains showed decreased COS -induced permeabil- 
ity compared to wild type cells. Because COS is a cell 
stress causing agent, we compared the transcriptional 
response induced by COS to those of other environ- 
mental stressors and ascertained whether an environ- 
mental stress pretreatment could provide resistance to 
COS. Finally, we evaluated whether a COS-resistant 
ARL1 overexpressing strain is also resistant to other 
antifungal compounds. 

Results 

Sensitivity to chitosan and COS-5.44 treatment in the 
yeast deletion pools 

A scheme summarizing the experiments and analyses 
performed in this study is shown in Figure 1. We first 
compared yeast sensitivity to deacetylated chitosan (ap- 
proximately 70 kDa, 80% deacetylated) and to chitosan 
oligosaccharide (COS) (5.44 kDa, 97% deacetylated). 
COS (Additional file 1: Figure SI A) had considerably 
higher antifungal activity than chitosan at equivalent 
concentrations (Additional file 1: Figure SIB). Because 
of this, it is easier to more consistently manufacture this 
form, we decided to use COS as the antifungal agent. 
The heterozygous and homozygous yeast deletion collec- 
tions were grown in YPD with 91.1 Hg/ml COS-5.44, a 
concentration that inhibited yeast wild type growth by 
10 - 20% compared with the control (Figure 2B). In the 
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Experimental approach 
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(See figure on previous page.) 

Figure 1 Experimental approach used to identify and characterize genes that confer resistance to chitosan oligosaccharide (COS). 

Three chemogenomic assays were used: Haploinsufficiency profiling (HIP), homozygous profiling (HOP) and multicopy suppression profile (MSP). 

(1) Heterozygous, homozygous deletion pools and multicopy suppression pool were grown competitively in the presence of COS-5.44. If a gene 

is required to grow in the presence of COS, the corresponding deletion strain will grow more slowly and therefore will be underrepresented. 

Cells overexpressing a gene that suppresses sensitivity will growth faster and will be overrepresented in the MSP pool. (2) Genomic DNA was 

isolated from cells prior to and after the HIP-HOP assays, and plasmid purification from the COS treated MSP pool was carried out. (3) Barcodes 

were PCR amplified for HIP-HOP assays as well as the plasmid inserts of MSP. (4) PCR products (barcodes and plasmid inserts) were hybridized to 

a TAG4 array. Intensity of treatment samples is compared with intensity of a control sample to determine relative abundance (~ fitness). (5) 

Sensitive deletion strains and constructed overexpressing strains were individually confirmed. (6) Five resistant overexpressing strains that were 

also sensitive as deletion strains were selected for transcriptome analysis. Overexpressing strains and vector control were grown in the presence 

of COS (1 1 2.5 ug/ml) and cells were harvested before COS treatment and after 60 min of COS treatment. (7) RNA was isolated from harvested 

cells, cDNA synthesized and labelled with fluorescent dyes. (8) Labelled samples were hybridized to expression microarrays. Transcriptional 

changes were indentified by differential expression analysis. Figure modified from Ericson et al. (2010). 
* ) 



HIP-HOP assays, 39 yeast deletion strains sensitive to 
COS-5.44 were indentified (22 homozygous and 17 het- 
erozygous; log2 ratio > 3.5, see Methods; Table 1, Add- 
itional file 2: Table SI). 

A previous global fitness analysis, similar to the HIP- 
HOP assay we performed, identified 101 chitosan sensi- 
tive homozygous and three heterozygous deletion 
strains, [23]. Approximately 10% of the homozygous de- 
letion strains found in their study were among the genes 
identified in our screen (ARL1, IMP2', APN1, RSC1, 
SNF8, DFG16, VPS66, YAP3, SSN8), but none of the 
three heterozygous deletions strains corresponded to 
those identified in our screen. 

Biological processes (defined by Gene Ontological 
terms) associated with the COS-5.44 sensitive deletion 
strains were determined using gene set enrichment ana- 
lysis (GSEA) [29]. The enriched processes included 
RNA biology (transcription, translation and regulatory 
mechanisms), membrane functions (e.g. transport and 
targeting), membrane structural components (e.g. pro- 
teins), cell division (spindle body and microtubules) and 
proteasome processes (structural and regulatory func- 
tions; Figure 2C, FDR < 0.1, see Methods). Similar GO 
terms showed altered fitness in the presence of chitosan 
in a similar fitness analysis [23]. 

To confirm whether the 39 deletion strains were sensi- 
tive to COS on their own, we tested each strain individu- 
ally and found that 21 (-50%) were inhibited by COS 
(Table 1). 

Resistance to COS-5.44 treatment in the multicopy 
suppression profiling pool 

MSP assays involve the transformation of yeast with a 
multicopy suppression plasmid library containing ran- 
dom fragments of approximately 7 kbp pieces of DNA. 
The MSP assay was performed with 101.25 - 250.0 ug/ 
ml COS-5.44 and it was determined that 112.5 ug/ml 
COS-5.44 inhibited wild type yeast growth by at least 
70% (Figure 2B). Yeast containing suppression plasmids 
that provide resistance to 112.5 ug/ml COS-5.44 will 



expand more quickly than the rest of the yeast of the 
population. Genomic DNAs from a resistant population 
and a transformed but untreated population are 
extracted and hybridized to ORF microarrays to identify 
genes present in multicopies in the resistant population. 

Two independent replicates of the MSP screen identi- 
fied a total of 68 genes as putative suppressors of COS- 
5.44 sensitivity (Additional file 1: Figure S2). The 
random genomic DNA fragments in the multicopy sup- 
pression plasmid library contain on average 2-3 yeast 
genes per fragment [25]. Only one gene in each frag- 
ment is likely responsible for resistance; therefore we 
expected only 33-50% of these genes to be confirmed as 
suppressors of COS-5.44 (see below and Additional file 
2: Table S2). There were 42 and 32 suppressing genes in 
replicates 1 and 2 respectively with only 6 genes in com- 
mon (Additional file 1: Figure S2). The 6 genes were: 
PKR1, a V-ATPase assembly factor in the endoplasmic 
reticulum (ER) which functions with other V-ATPase as- 
sembly factors to efficiently assemble the V-ATPase 
membrane sector; COX5A, subunit Va of cytochrome c 
oxidase, involved in mitochondrial electron transport; 
FLOl, a lectin-like cell wall protein involved in floccula- 
tion that binds to mannose on the surface of other cells; 
VAC7, a vacuolar membrane protein involved in inherit- 
ance and vacuole morphology; RPL20A, a large riboso- 
mal subunit protein; and MSG5, a dual specificity 
phosphatase protein required for maintenance of low 
level signaling through the cell integrity pathway which 
regulates and is regulated by Slt2p. Three of these genes, 
PKR1, VAC7 and MSGS, were also found to be sensitive 
to COS when tested as single heterozygous deletion 
strains. Among the suppressing genes found in replicate 
2, ARL1, a gene encoding a GTPase protein, was identi- 
fied. This gene was found to be sensitive to COS-5.44 as 
a deletion strain in the HIP-HOP assay (Table 1). 

Confirmation of COS suppressing strains 

After the identification of putative resistant strains we 
tested each strain individually. From the 68 candidate 
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Figure 2 Chitosan oligosaccharide molecule, screen and biological process associated with sensitive deletions strains. A) Chitosan 
molecular structure. B) Chitosan oligosaccharide screen in wild type yeast (BY4743) at 5 concentrations (91.1 - 1 12.5 ug/ml) in 0.5X YPD pH 5. 
Optical density readings were taken every 15 min over 20 hrs using a Tecan Genios reader (see Methods). C) Biological processes associated with 
chitosan oligosaccharide (COS) sensitive deletion strains. A node represents a biological process significantly enriched in the COS-5.44 HIP-HOP 
assay (FDR < 0.1, see Methods). The node size correlates to the number of genes annotated to that functional category. The width of the edge 
correlates to the degree of gene overlap between the 2 connected categories. If the overlap coefficient is less than 0.5, edges are not shown 
(see Methods). Cluster membership is shown by node color, where clustering is based in degree of overlap among categories. Bar plots show the 
sensitivity scores of the genes that contributed to the functional enrichment of the cluster, border color surrounding the plot correlated with the 
nodes in the cluster. The top 10 genes in each cluster category are shown in each plot. Heterozygous strains are marked with a black dot. 
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Table 1 Top 39 characterized sensitive deletion strains found in the HIP-HOP assays above 3.5 log 2 ratio 



probeid 



COS-5.44 
(91.1 ug/ml) 
(log2ratio) 



Gene 



SGD_description 



Fitness 
score (avg) 2 



Deletion 
pool source 3 



YBR 164C 



YBR257W 



YIL048W 
YIL075C 
YIL154C 

YJL002C 



3.534 



YDL060W 3.769 

YDL213C 5.101 

YDR300C 4.371 

YDR320C 3.571 

YDR324C 3.636 

YDR389W 4.281 



YDR455C 3.577 
YGL045W 4.229 



YGR056W 3.72 



YGR122W 3.719 



YGR246C 3.887 



YHL009C 3.625 
YIL019C 3.59 

YIL022W 3.931 



3.83 

3.653 

4.315 

3.515 



ARL1 Soluble GTPase of the RAS superfamily, with a role in 0.756 HOM 

regulation of membrane traffic mainly regulates 
pottassium influx; similar to ADP-ribosylation factor. 

P0P4 subunit of both RNase MRP, which cleaves pre-rRNA, 0.9963 HET 
and nuclear RNase P, which cleaves tRNA precursors 
to generate mature 5'ends; binds to the RPR1 RNA 
subunit in RNase P. 

TSR1 Protein required for processing of 20S pre-rRNA in the 0.8887 HET 
cytoplasm associates with pre-40S ribosomal particles. 

NOP6 rRNA-binding protein required for 40S ribosomal subunit 0.8854 HOM 

biogenesis; contains an RNA recognition motif (RRM). 

PR01 Gamma-glutamyl kinase, catalyzes the first step in proline 0.8521 HOM 

biosynthesis. 

SWA2 Auxilin-like protein involved in vesicular transfort; clathrin- 0.9621 HOM 

binding protein required for uncoating of clathrin- 
coated vesicles 

UTP4 Subunit of U3 involved in production of 18S rRNA, 0.761 HET 

assembly of small ribosomal subunit and transcription 
of 35S rRNA transcript, member of t-Utp subcomplex. 

SAC7 GTPase activating protein (GAP) for Rholp, involved 0.8768 HOM 

in signaling to the actin cytoskeleton, null mutations 
suppress tor2 mutations and temperature sensitive 
mutations in actin; potential Gdc28p substrate. 

YDR455C Partially overlaps the verified gene YDR456W. 0.9312 HOM 

RIM8 Protein involved in proteolytic activation of Rim 101 p 0.896 HOM 

in response to alkaline pH; interacts with ESCRT-1 subunits 
Stp22p and Vps28p; member of the arrestin-related 
trafficking adaptor family. 

RSG1 Component of the RSC chromatin remodeling complex; 0.8656 HOM 

required for expression of mid-late sporulation-specific 
genes. 

YGR122W Probable ortholog of A. nidulans PalC, which is involved 0.9913 HOM 

in pH regulation and binds to the ESCRT-lll complex; 
null mutant does not properly process Rim 1 01 p and 
has decreased resistance to rapamycin. 

BRF1 TFFIIIB B-related factor, one of three subunits of RNA 0.9782 HET 

polymerase III transcription initiation factor TFFIIIB, 
binds TFIIIC and TBP and recruits RNA pol III to promoters, 
amino-terminal half is homologous to TFIIB. 

YAP3 Basic leucine zipper (bZIP) transcription factor. 0.8854 

FAF1 Protein required to pre-rRNA processing and 40S 0.9603 HET 

ribosomal subunit assembly. 

TIM44 Essential component of the Translocase of the inner 0.8167 HET 

Metochondrial membrane; tethers the import motor 
and regulatory factors (PAM complex) to the translocation 
channel (Tim23p-Tim1 7p core complex). 

NE01 Putative aminophospolipid translocase (flippase) involved 1.0734 HET 

in endocytosis and vacuolar biogenesis. 

RPN2 Subunit of the 26Sproteasome, substrate of the 0.9832 HET 

N-acetyltransferase Natlp. 

IMP2 Transcriptional activator involved in maintenance of 0.2516 HOM 

ion homeostasis and protection against DNA damage 
caused by bleomycin and other oxidants. 

OST1 Alpha subunit of the oligosaccharyltransferase 1.0555 HET 

complex of the ER lumen, which catalyzes asparagine- 
linked glycosylation of newly synthesized proteins. 
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Table 1 Top 39 characterized sensitive deletion strains found in the HIP-HOP assays above 3.5 log 2 ratio (Continued) 



YGL188C 
YJR102C 
YKL114C 

YKR024C 
YKR025W 

YLR100W 

YLR147C 
YLR223C 

YLR378C 

YNL025C 

YNL287W 
YNR046W 

YOR014W 
YOR030W 

YOR080W 

Y0R117W 

YOR244W 



3.656 



3.755 



4.271 



3.506 
4.622 



3.929 



3.525 



3.633 



3.539 



3.581 



4.047 



3.844 
3.669 



3.506 



3.5 



BUD19 ORF overlaps the verified gene RPL39 by 88%f; diploid 1.0087 HOM 

mutants displays a weak budding pattern phenotype 
in a systematic assay. 

VPS25 Component of the ESCRT-II complex, which is involved 1.0011 HOM 

in ubiquitin -dependent sorting of proteins into the 
endosome. 

APN1 Apurinic/apyrimidinic endonuclease, 3'-repair diesterase 0.8878 HOM 

involved in repair of DNA damage by oxidation 
and alkylating agents; also functions as a 3-5' exonuclease 
to repair 7,8-dihydro-8-oxodeoxyguanosine. 

DPB7 Putative ATP-dependent RNA helicase of the dead 0.9795 HOM 

-box family involved in ribosomal biogenesis. 

SNF7 One of four subunits of the endosomal sorting complex 0.8599 HOM 

required for transport III (ESCRT-III); involved in the 
sorting of transmembrane proteins into the multivesicular 
body (MVB) pathway. 

ERG27 3-keto sterol reductase, catalyzes of the last three steps 0.8809 HET 

required to removed two C-4 methyl groups from an 
intermediate in engosterol biosynthesis; mutants are 
sterol auxotrophs. 

SMD3 Core Sm protein Sm D3; part of heteroheptameric complex 1.0582 HET 

that is part of the spliceosomal U1, U2, 114, and U5 snRNPs; 
homolog of human Sm D3. 

IFM1 Coactivator that regulates transcription of ribosomal protein 0.9691 HET 

(RP) genes; recruited to RP gene promoters during optimal 
growth conditions via Fhllp; subunit of CURI, a complex 
that coordinates RP production and pre-rRNA processing. 

SEC61 Essential subunit of Sec61 complex (Sec61p, Sbhlp, and 0.8717 HET 

Ssslp); with Sec63 complex forms a channel for SRP-dependent 
protein import and retrograde transport of misfolded 
proteins out of the ER. 

SSN8 Cyclin-like component of the RNA polymerase II holoenzyme, 0.823 HOM 

involved in phosphorylation of the RNA polymerase II 
C-terminal domain, in glucose repression and telomere 
maintenance. 

SEC21 Gamma subunit of coatomer, a heptameric protein complex 0.8959 HET 

that together with Arflp forms the COPI coat; involved in ER 
to Golgi transport of selective cargo. 

TRM112 Subunit of tRNA methyltransferase (MTase) complexes in 0.9906 HET 
combination with Trm9p and Trm11p; subunit of complex 
with Mtq2p that methylates Sup45p (eRFI) in the ternary 
complex eRF1-eRF3-GTP. 

RTS1 B-type regulatory subunit of protein phospatase 2A (PP2A); 0.9832 HOM 

homolog of the mammalian B'subunit of PP2A. 

DFG16 Probable multiple transmembrane protein, involved in 1.0307 HOM 

diploid invasive and pseudohyphal growth upon nitrogen 
starvation; required for accumulation of processed RimlOlp. 

DIA2 Origin-binding F-box protein that forms an SCF ubiquitin 1.1312 HOM 

ligase complex with Skpl p and Cdc53p; plays a role in DNA 
replication, involved in invasive and pseudohyphal growth. 

RPT5 One of the six ATPases of the 19S regulatory particle 0.8955 HET 

(26S proteasome) involved in the degradation of 
ubiquintinated substrates; recruited to the GAL1-10 promoter 
region upon induction of transcription; similar to human TBP1. 

ESA1 Catalytic subunit of the histone acetyltransferase complex 0.8716 HET 

(NuA4) that acetylates four conserved internal lysines of histone 
H4N-terminal tail; required for cell cycle progression and 
transcriptional silencing at the rDNA locus. 
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Table 1 Top 39 characterized sensitive deletion strains found in the HIP-HOP assays above 3.5 log 2 ratio (Continued) 



YPL002C 


3.679 


SNF8 


Component of the ESCRT-II complex, which is involved in 
ubiquitin-dependent sorting of proteins into the endosome; 
appears to be functionally related to SNF7; involved in 
glucose derepression. 


0.8925 


HOM 


YPR139C 


3.643 


VPS66 


Cytoplasmic protein of unknown function involved in 
vacuolar protein sorting. 


0.8768 


HOM 



1 COS-5.44 log2 ratio is calculated from tag4 array data for each tag of each deletion strain. COS-5.44 log2 ratio = log2 (u c - bg)/([i t - bg), where u c is the mean 
intensity of the control sample, u t is the mean intensity of the treated samples, and bg is the mean intensity of the unassigned probes. This score is proportional 
to the log2 ratio of cells present in the control sample versus the treatment sample (i.e. deletion strains exposed to 91.1 ug/ml of COS) [24]. 

2 Fitness score (W) is estimated from the average growth curves (as measured by optical density) of a given strain in the presence and absence of COS after 20 h. 
The fitness scores from three individual colonies of each deletion strain grown in triplicate (i.e. 9 cultures for each strain) were used to calculate the average. This 
score takes into account the growth of the wild type and deletion strain in control conditions as well as treatment conditions. lV c = wt c /st c ; Wt = wt t /st t , where W c 
is the fitness score under control (no drug) conditions, W t is the fitness score after treatment, wt is doubling time of wild type strain, st is doubling time of the 
deletion strain, c is control conditions, t is treatment (COS) conditions. Normalized fitness (tV) = avg W t /avg W c [24]. A fitness score less than 1 indicates that the 
strain is sensitive [24]. A minimum difference of 0.1 (10%) was used as a cut-off to consider a strain sensitive. 

3 Deletion pool source identifies which pool of deletion strains the deleted strain came from. HET - heterozygous essential deletions strains, HOM - homozygous 
deletion strains. 



suppressor genes, 57 individual overexpressing/suppressing 
strains each containing a single putative resistance gene 
were obtained from a collection of ORFs whose expres- 
sion is driven by native promoters [30] see Methods. 
Twenty-one (-31%) of the putative suppressor genes 
conferred resistance to COS-5.44 when tested in this 
manner (Table 2). COS-5.44 resistant overexpressing 
strains were able to grow in 112.5 ug/ml COS-5.44 while 
the wild type (transformed with an empty vector) was un- 
able to grow (Figure 3B). 

Among the confirmed genes that confer resistance to 
COS-5.44 were two essential genes: RBAS0 (YDR527W), 
a gene whose protein product is involved in transcrip- 
tion, and ARP2 (YDL029W), a gene encoding actin 
related protein 2, a subunit of the Arp2/3 complex 
which is required for the motility and integrity of cor- 
tical actin patches. Among the remaining 19 non- 
essential genes, a diversity of functions were found, such 
as a sterol reductase, a heat shock protein, ADP- 
ribosylation factors, soluble GTPases, alpha tubulin, and 
a ubiquitin conjugation enzyme. Putative suppressors of 
COS antifungal action in yeast were all tested as hetero- 
zygous deletion strains and from the 21 resistant overex- 
pressing strains, 11 were found to be sensitive when 
tested as heterozygous deletion strains (Table 2). Five 
genes that provided COS resistance when overexpressed 
or sensitivity as a heterozygous or homozygous deletion 
strain were selected for further study primarily based 
on the known or putative functions of the genes. More 
specifically we selected genes with roles in signalling 
pathways, cell or membrane integrity, or transcription 
regulation. 

We selected the following 5 overexpressing strains: 
ARL1, which encodes a GTPase involved in membrane 
trafficking, was selected because it was found in the 
HIP-HOP assay as a sensitive homozygous deletion 
strain and in the MSP assay as a multicopy suppressor. 
The other overexpressing strains we selected were: 



BCK2 and MSG5, which are both involved in cell integ- 
rity pathways; ERG24, a gene involved in ergosterol syn- 
thesis; and RBAS0 (mentioned above). BCK2 is a Ser- 
Thr rich protein with protein kinase C activity that acts 
in signal transduction. Overexpression of BCK2 can res- 
cue defects in a cwh43A mutant that displays several cell 
wall defects [31]. BCK2 overexpression can also suppress 
the cell lysis defects seen when the kinases Mpkl and 
Pckl are deleted [32]. MSG5 encodes a protein phos- 
phatase involved in cell cycle control through the depho- 
sphorylation of MAPK and is required for restricting 
signaling by the cell integrity pathway in yeast [33]. The 
inhibition of MAPK signaling leads to inhibition of cell 
differentiation and cell division [34]. The functions of 
ARL1 and ERG24 and their potential roles in chitosan 
resistance are described in more detail in the Discussion. 

The fitness of all overexpressing strains compared to 
an empty vector control and the corresponding hetero- 
zygous deletion strains were assessed. In the presence of 
COS-5.44, the overexpressing strains always grew better 
than the vector control (Figure 3B, showing results of 
ARL1) and the corresponding deletion strain demon- 
strated reduced fitness (Figure 3D). In the absence of 
the COS-5.44, no growth differences were observed 
among each overexpressing strain or deletion strains 
(tested in triplicate) compared with the vector control, 
indicating the absence of inherent fitness defects 
(Figure 3A and 3C). 

Transcriptional response to COS 

To identify the transcriptional changes caused by expos- 
ure to COS-5.44, wild type (vector control) cells, were 
exposed to COS for 1 h, total RNA extracted, and a 
microarray analysis was performed. A total of 335 genes 
were differentially expressed in response to COS-5.44 
(P-value < 0.05 and log2 fold change > 1 or < -1). 

A GSEA was performed for the entire transcription 
data set and among the up-regulated transcripts, we 



Table 2 Twenty-one yeast overexpression strains confirmed as suppresors of sensitivity to COS-5.44 



ORF 



Gene Resistance Fitness score of 
confirmed resistance (avg) 3 



Sensitivity confirmed 
(HET deletion strain) 



Fitness score of 
sensitivity (avg) 4 



Description from SGD 



YBR164C 

YBR166C 
YDL029W 

YDR171W 



YLR285W 
YML124C 
YMR123W 
YNL053W 



ARL1 1 

TYR1 

ARP2 

HSP42 



YDR524C AGE1 

YDR527W RBA50 1 

YER048C CAJ1 

YER167W BCK2 1 

YJL046W AIM22 

YKL208W CBT1 

YLR193C UPS1 



NNT1 
TUB3 
PKR1 2 
MSG5 1 



YNL054W VAC7 2 



yes 

yes 
yes 

yes 

yes 
yes 
yes 

yes 

yes 
yes 
yes 

yes 
yes 
yes 
yes 

yes 



1.870 

1.381 
1.609 

1.372 

1.783 
1.407 
1.845 

1.365 

1.258 
1.265 
1.201 

1.448 
1.763 
1.350 
1.454 

1.387 



yes 



yes 



yes 



yes 
yes 

yes 



0.756 

1.035 
1.023 

1.135 

1.106 
0.899 
0.995 

0.890 

0.908 
0.930 
0.953 



1.227 
0.866 
0.874 

0.873 



Soluble GTPase of the Ras superfamily, with a role in regulation of 
membrane traffic mainly regulates potassium influx; 
similar to ADP-ribosylation factor. 

Prephenate dehydrogenase involved in tyrosine biosynthesis, 
expression is dependent on phenylalanine levels. 

Essential component of the Arp2/3 complex involved in endocytosis 
and membrane growth and polarity. A conserved actin nucleation 
center required for the motility and integrity of actin patches. 

Small heat shock protein (sHSP) with chaperone activity involved 
in cytoskeleton reorganization after heat shock; forms barrel-shaped 
oligomers that suppress unfolded protein aggregation. 

ADP-ribosylation factor (ARF) GTPase activating protein (GAP) effector, 
involved in the secretory and endocytic pathways. 

Protein involved in transcription; interacts with RNA polymerase I 
subunits Rpb2p, Rpb3, and Rpb1 1 p; has similarity to human RPAP1 . 

Nuclear type II J heat shock protein of the of the E. coli dnaJ family, 
binds, to non-native substrates for presentation to Ssa3p, may function 
during protein translocation, assembly and disassembly. 

Protein rich in serine and threonine residues involved in protein kinase 
C signaling pathway, which controls cell integrity; overproduction 
suppresses pkc 1 mutations 

Putative lipoate-protein ligase, required along with Lip2 and Lip5 for 
lipoylation of Latlp and Kgd2p. 

a role in 3'end processing of the COB pre-mRNA; displays genetic 
interaction with cell cycle-regulated kinase. 

Mitochondrial intermembrane space protein that regulates mitochondrial 
cardiolipin levels, null has defects in Mgmlp processing, integrity of 
mitochondrial inner membrane complexes; ortholog of human PRELI. 

Putative nicotinamide N-methyltransferase, has a role in rDNA silencing 
and in the lifespan determination. 

Alpha-tubulin; associates with beta-tubulin (Tub2p) to form tubulin dimer, 
which polymerizes to form microtubules; expressed at lower level than Tub1p 

V-ATPase assembly factor, functions with other V-ATPase assembly factors 
in the ER to efficiently assemble the VATPase membrane sector (VO). 

Dual-specificity protein phosphatase; required for maintenance of a low 
level of signaling through the cell integrity pathway, adaptive response to 
pheromone; regulates and is regulated by Slt2p; dephosphorylates Fus3p. 

Integral vacuolar membrane protein involved in vacuole inheritance and 
morphology; activities Fablp kinase activity under basal conditions anfd 
also after hyperosmotic shock. 



Table 2 Twenty-one yeast overexpression strains confirmed as suppresors of sensitivity to COS-5.44 (Continued) 



YNL218W 


MSG1 2 


yes 


1.214 


yes 


0.899 


Protein with DNA-dependent ATPase and ssDNA annealing activities 
involved in maintenance of genome; interacts functionally with DNA 
polymerase delta; homolog of human WHIP. 


YNL280C 


ERG24 1 


yes 


1.275 


yes 


0.893 


C-14 sterol reductase, acts in ergosterol biosynthesis; mutants accumulate 
the abnormal sterol ignosterol (ergosta-8, 14dienol). 


YNR057C 


BI04 2 


yes 


1.246 


yes 


0.885 


Dethiobiotin synthetase, catalyzes the third step in the biotin biosynthesis 
pathway; BI04 is in a cluster of 3genes (BI03, BI04, and BI05) that mediate 
biotin synthesis. 


YPL053C 


KTR6 2 


yes 


1.616 


yes 


0.861 


Probable mannosylphosphate transferase involved in the synthesis of core 
oligosaccharides in protein glycosylation pathway; member of the 
KRE2/MNT1 mannosyltransferase family. 


YPL106C 


SSE1 2 


yes 


1.253 


yes 


0.896 


ATPase that is a component of the heat shock protein Hsp90 chaperone 
complex; binds unfolded proteins; member of the heat shock protein 70 
(HSP70) family; localized to the cytoplasm. 



1 Genes selected for transcriptomic analysis. 

2 Overexpressing strains confirmed as suppressors when overexpressed and sensitive when tested as deletion strains. 

3-4 Fitness score (W) is estimated from the average growth curves (as measured by optical density) of a given strain in the presence and absence of COS after 20 h. Data shown corresponds to deletion strains exposed 
to 91.1 ug/ml and overexpressing strains to 1 12.5 ug/ml of chitosan oligosaccharide. The fitness scores from three individual colonies of each overexpressing or deletion strain grown in triplicate (i.e. 9 cultures for 
each strain) were used to calculate the average. This score takes into account the growth of the wild type and overexpressing or deletion strain in control conditions as well as treatment conditions. W c = wt c /st c ; 
Wt = wt t /st t , where W c is the fitness score under control (no drug) conditions, W t is the fitness score after treatment, wt is doubling time of wild type strain, st is doubling time of the overexpressing or deletion strain, c 
is control conditions, t is treatment (COS) conditions. Normalized fitness (W) = avg W t /avg W c [24]. A fitness score less than 1 indicates that the strain is sensitive, a fitness score above 1 indicates that the strain is 
resistant [24]. Means that the strain has a treatment induced growth defect, a minimum difference of 0.1 (10%) was used as a cut-off to consider a strain sensitive. 
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Figure 3 Confirmation of COS-5.44 resistance in the Arl1 overexpressing strain and sensitivity of the corresponding heterozygous 
deletion strain. The overexpression of ARLl does confer resistance to COS-5.44 (1 12.5 ug/ml) compared with the wild type (vector control) that 
is not able to grow at this concentration. The heterozygous deletions strain (aril A) shows sensitivity to COS-5.44 (91.1 ug/ml). A) Growth curves 
of ARLl and wild type (vector control) growth in the vehicle. B) Growth curves of Aril overexpressing and wild type strains growth with COS-5.44 
(1 12.5 ug/ml). C) Growth curves of the heterozygous deletion (aril A) and wild type strain grown with the vehicle. D) Growth curves of 
heterozygous deletion (aril A) and wild type strains grown with COS-5.44 (91.1 ug/ml). Optical density readings were taken every 15 min over 
20 hrs using a Tecan Genios reader. Tecan ODs were converted to conventional 1 mm path length cuvette ODs using a calibration function 
provided by Ericson et al. 201 0, [24]. Similar behaviour was observed for the rest of the selected overexpressing strains (Bck2, Erg24, Msg5 and 
Rba50) deduced to be resistant to chitosan oligosaccharide COS-5.44 from the MSP screen. The heterozygous deletion strains from these genes 
(bck2A, erg24A msg5A and rba50A) were also individually tested for sensitivity to COS-5.44 (91.1 ug/ml). Three colonies of each strain were 
grown in triplicate and compared with the wild type in the presence of the COS-5.44 and the vehicle (1% DMSO). 



identified cellular respiration, ATP production, protein 
complex biogenesis, and mitochondria translation and 
organization as enriched biological processes (P< 0.005). 
For the down-regulated transcripts, the enriched terms 
included glycosylation, transmembrane transport, and 
sterol and lipid biosynthesis. For the overexpressing 
strains that provided resistance to COS (see below), the 
majority of these biological processes had the opposite 
characteristic (i.e. enriched biological processes among 
up-regulated transcripts in the wild type were enriched 
among down-regulated transcripts in the overexpression 
strains). A comparison of the transcriptional profiles of 
the COS-5.44 treated wild type cells in this study with 
the chitosan treated 60 min profile from another similar 
study, [22] was performed (Additional file 1: Figure S3, 
Additional file 2: Table S13). As might be expected, there 
was considerable overlap among enriched biological pro- 
cesses between the two studies (Additional file 1: Figure 
S3, Additional file 2: Table S13), and included such 
terms as cell wall organization, ATP production, and 



oxidative phosphorylation as being enriched among up- 
regulated transcripts. Ribosome biogenesis and poly- 
phosphate metabolic processes were enriched among 
down-regulated transcripts. There were also several bio- 
logical processes that showed up as having opposite 
enrichments in the two studies. For example, sterol/ 
ergosterol/lipid biosynthesis was observed as a down- 
regulated process in our study while these were up- 
regulated in Zakrzewska and co-workers study. While it 
is not clear why some processes show opposite enrich- 
ments in the two studies, we do know that different 
types and amounts of anti-fungal agents were used, 
Zakrzewska et al. exposed the yeast cells to fragmented 
chitosan at an ICi 0 . 15 and in this study, cells were 
exposed to COS-5.44 at an IC 70 _ 80- 

Transcriptome analysis of the COS-5.44 resistant strains 

To gain a further understanding on the mode of action 
and mechanisms of resistance to COS, we performed a 
transcriptome analysis of the 5 overexpressing strains 
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Figure 4 (See legend on next page.) 
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(See figure on previous page.) 

Figure 4 Differentially expressed genes in COS-5.44 resistant yeast overexpressing strains in the presence of COS-5.44. One thousand 
two hundred and twenty genes showed a significant change in expression in at least one of the overexpressing strains compared with the wild 
type (P-value<0.05 and log2 fold change > 1 or < — 1 ). Clustering of the 1220 genes with significant change in expression was based in similarity 
(see methods). A-D) Subsets of up-regulated genes among the 5 overexpressing strains. E-H) Subsets of down-regulated genes among the 5 
overexpressing strains. 



known to increase resistance to COS-5.44 (data available 
in GEO under GSE32888). Each overexpressing strain 
and the wild type (vector control) were treated with 
COS-5.44 for 1 h or mock treated with vehicle alone 
and RNA was isolated from each sample. The RNA was 
converted to labeled cDNA and hybridized to Nimble- 
Gen expression microarrays (see Methods). We per- 
formed two sets of analyses. The first analysis was 
designed to find genes differentially expressed genes in 
each overexpression strain compared to wild type in the 
absence of COS in order to see if overexpression of the 
gene in question changed the transcriptional profile (and 
potentially the physiology or structure of the cell). The 
second analysis was designed to find differentially 
expressed genes in each overexpression strain compared 
to wild type in the presence of COS to see if the over- 
expression strains responded differently to COS than 
wild type. 

In the first analysis, we identified 184 genes with differ- 
ential expression in at least one of the 5 overexpressing 
strains compared with the wild type (vector control) with- 
out COS-5.44 treatment (Additional file 1: Figure S4, 
Additional file 2: Table S3; P-value<0.05 and log2 fold 
change >1 or<-l). Each of the 5 overexpressing strains 
showed increased expression of the gene contained in the 
transformation vector (Additional file 1: Figure S4 A-D). 
There was a subset of 13 highly up-regulated genes in the 
ARL1 and RBA50 strains (Additional file 1: Figure S4 A) 
including genes encoding heat shock proteins, acid phos- 
phatases, inorganic phosphatases, and transmembrane 
transporters. Several genes involved in oxidative phos- 
phorylation, amino acid and carbohydrate metabolic pro- 
cesses and biosynthesis, and mitochondrial ATP synthesis 
and electron transport were up-regulated (Additional file 
1: Figure S4 D). 

Transcriptome analysis following COS-5.44 treatment 
was performed and 1220 genes (-589 up- and -631 
down-regulated) showed significant differential expres- 
sion relative to the wild type under the same growth 
conditions (P-value < 0.05 and log2 fold change > 1 or 
< -1, Additional file 2: Table S4) in at least one of the 5 
overexpressing strains (Figure 4). Significantly over 
represented biological processes identified by enrich- 
ment analyses [35,36] in the 589 differentially up- 
regulated genes included transcription, cell cycle, protein 
modification, stress response and RAS signal transduc- 
tion (Figure 4A-E). Down-regulated genes were enriched 



for biological processes included protein folding, protein 
complex assembly, and respiratory chain complex genes 
(Figure 4F-H). 

The primary biological processes associated with the 5 
overexpression strains in the presence of COS (Table 2) 
were: membrane signalling functions (ARL1, BCK2, 
MSGS), transcription {RBA50) and ergosterol synthesis 
(ERG24). To gain insight into the potential mechanisms 
of COS resistance for each overexpressed gene, enrich- 
ment maps were constructed from the entire transcrip- 
tion data set, GSEA analyzed, for each of the COS 
treated overexpressed strains (Figure 5 and Additional 
file 1: Figures S5- S9). Genes that were down-regulated 
to various degrees in all 5 strains were associated with 
processes such as cell energy generation (mitochondrial 
biology, ATP metabolism, energy storage metabolites) 
and associated by-products (oxidative stress). Most over- 
expressing strains (ARL1, ERG24 and RBA50) displayed 
up-regulated genes involved in cell cycle progression 
(mitosis/meiosis, chromatin dynamics and modification 
and sporulation) and transcription. Taken together 
these results suggest that the overexpression resistant 
strains had an overall reduction in energy production 
and an increase in cell proliferation in response to 
COS-5.44 perturbation compared to wild type cells. 

We confirmed some of the global transcriptional 
changes by qRT-PCR. ARL1 overexpression was con- 
firmed in Aril overexpressing strains. For all of the 
genes selected, the qRT-PCR results were always in the 
same direction as the microarray results (i.e. increased 
transcript levels for a given gene found using microar- 
rays was also found to be increased by qRT-PCR) but 
the dynamic range for the transcript changes was 
greater in the qRT-PCR assays. COS-5.44 treatment 
increased levels of expression in the Aril overexpres- 
sing strain for FIG2, MUC1, VBA5 and YJU2 with the 
qRT-PCR transcript levels being 50 - 200% higher 
when compared with the microarray data. Decreased 
transcript levels were confirmed for CMD1, COX5B, 
HSP30, RCR1 and UBI4, using qRT-PCR (Additional 
file 1: Figure S10). 



ARL1 overexpression reduces COS-induced 
membrane permeabiiization 

A Sytox permeability assay was performed in the wild 
type (BY4743) and ARL1 strain after overnight growth in 
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Figure 5 Biological processes associated with differentially expressed genes on Aril overexpressing strain when exposed to chitosan 
oligosaccharide (COS). A node represents a biological process significantly enriched (FDR < 0.1, see Methods). Boxes on the side show summary 
of the main biological process found in a cluster. The node size correlates to the number of genes annotates to that functional category. Red and 
blue node border colours indicate enrichment of biological processes in up- and down-regulated genes, respectively. The width of the edge 
correlates to the degree of gene overlap between the 2 connected categories. If the overlap coefficient is less than 0.5, edges are not shown (see 
Methods). Cluster membership is shown by node color, where clustering is based in degree of overlap among categories. 
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Figure 6 Sytox cell permeability assay of COS-5.44 resistant 
Aril overexpressing strain vs. vector control (BY4743). A) Vector 
control (BY4743) shows significantly higher cell membrane 
permeability compared with overexpressing strain Aril (B) after 
treatment with chitosan oligosaccharide (COS-5.44). A large increase 
in cell membrane permeability is observed in the wild type as 
chitosan oligosaccharide concentration is increased (A). The first 
peak (left) in the graph corresponds to the background fluorescence 
and the second peak (right) corresponds to the Sytox signal 
observed inside the cell due to nucleic acid binding. C) Overlay of 3 
replicates of each strain Sytox assay after exposure to1 12.5 ug/ml 
COS-5.44. Student f-test of C) data, P-value < 0.01 01 . 



YPD with COS-5.44 to assess cell membrane permeabil- 
ity [37] see Methods. Sytox Green fluoresces once bound 
to nucleic acids. It only enters cells with compromised 
plasma membranes and is excluded from live cells with 
an intact plasma membrane. Wild type cells treated with 
COS-5.44 were Sytox-positive indicating that COS treat- 
ment causes membrane permeability (Figure 6A and 
6C). There was a dose dependent increase in the intra- 
cellular Sytox fluorescence signal as a function of COS 
concentration (Figure 6A). Although there was some 
Sytox signal in the Aril overexpressing strain, it was sig- 
nificantly lower than that in the wild type at the highest 
concentration of COS tested (112.5 ug/ml; P-value < 0.011, 
student t- test) (Figure 6C). This indicates that overex- 
pression of Aril provides protection against COS-induced 
cell permeability and damage. 



Yeast stress responses and COS-5.44 resistance 

Given our observation that the transcriptional changes 
in the overexpressing strains include genes involved in 
stress responses, we investigated whether these strains 
had activated well known stress responses that could po- 
tentially account for the strain's resistance to COS-5.44 
treatment. To test this hypothesis, the wild type strain 
(BY4743) was exposed to a sub-lethal dose of different 
primary stresses previously reported by Berry and Gasch 
2008, [38] followed by exposure to acute stress. 

The wild type strain was first exposed to sub-lethal 
doses of stress (mild primary stress) followed by 2 h ex- 
posure to COS-5.44 (acute secondary stress). The pri- 
mary mild stresses tested were: 0.25 M and 0.5 M NaCl, 
1 M sorbitol, 0.003% and 0.006% H 2 0 2 for 60 min, 37°C 
for 15 min and 30°C control (see Methods). 

Wild type cells did not acquire resistance to COS-5.44 
treatment after exposure to sub-lethal doses of different 
primary stressors (Additional file 1: Figure Sll), al- 
though there was some variation in how the different 
primary stressors affected yeast growth when challenged 
with COS-5.44. In the 30°C control treatments, there 
was a similar amount of growth in both the control 
medium (0.5X YPD pH 5) and vehicle (1% DMSO); 
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Figure 7 Effect of osmotic stressors, antifungal agents and oxidative stress on the growth of Aril overexpressing and wild type cells. 

A) Aril overexpression confers resistance to COS-5.44 at a concentration that inhibits growth of the vector control cells. The Aril overexpressing 
strain is as sensitive as the wild type cells to the tested antifungal compounds: B) Amphotericin B; C-D) Fluconazole; E-F) Terbinafine. 
Overexpression of ARL1 confers tolerance to H 2 0 2 (G) compared with the wild type. All assays were done in YPD except COS-5.44 that was done 
in 0.5X YPD. Optical density readings were taken every 15 min over 20 hrs using a Tecan Genios reader. Tecan ODs were converted to 
conventional 1 mm path length cuvette ODs using a calibration function provided by Ericson et al. 2010, [24]. Three colonies of each strain were 
grown in triplicate and compared with the wild type grown under the same conditions. YPD - Yeast Peptone Dextrose Broth. 



however, when COS-5.44 treatment was applied, we 
observed growth inhibition (Additional file 1: Figure 
SUA). After stressing the wild type cells for 15 min at 
37°C, they could withstand 91.1 ug/ml of COS-5.44 
and had a similar growth rate as wild type cells, but in- 
creasing COS concentrations dramatically inhibited 
growth rates (Additional file 1: Figure SI IB). When the 
highest concentration of NaCl (0.5 M) was used as the 
primary stress, the effect of 112.5 (-ig/ml of COS-5.44 
was increased, similar to what was observed when stres- 
sing cells for 15 min at 37°C (Additional file 1: 
Figure S11C). The addition of H 2 0 2 slowed growth of 



the wild type cells and H2O2 pre-treatment did not ap- 
pear to provide resistance to COS-5.44 (Additional file 1: 
Figure S11E-F). 

In contrast with the other stressors, sorbitol (1 M) 
stress decreased growth in both untreated and COS-5.44 
challenged cells (Additional file 1: Figure S11D). It 
should be noted that sorbitol primary stress does appear 
to provide some protection against the higher concen- 
trations of COS with cells growing as well in the higher 
as the lower concentrations of COS, a property not 
observed for the other stressors. The osmotic stress gen- 
erated by sorbitol treatment is known to activate other 
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Figure 8 COS-5.44 has a synergetic inhibitory effect on cell 
growth when used in combination with Fluconazole. A) Growth 
curves of wild type strain (BY4743) growth in the vehicle, COS-5.44 
(77 pg/ml), Fluconazole (16 ug/ml) and COS-5.44 + Fluconazole 
(77 + 16 ug/ml, respectively). B) Growth curves of wild type strain 
(BY4743) growth in the vehicle, COS-5.44 (84 ug/ml), Fluconazole 
(20 pg/ml) and COS-5.44 + Fluconazole (84 + 20 pg/ml, respectively). 
Optical density readings were taken as described in Figure 3. 



pathways such as MAPK and HOG cascades [39] and 
our results show that sorbitol pre-treatment provides 
some protection against COS (Additional file 1: Figure 
SI ID) [39]. The sorbitol results are consistent with the 
findings of Zakrzewska et al. 2007, [23] who found that 
when HOG pathway mutants and wild type cells, were 
exposed to 1 M sorbitol, partial protection against chito- 
san was observed. Overall, these results suggest that the 
COS-5.44 stress response is different from the previously 
described environmental stress responses and as a result, 
most pre-stresses do not confer resistance to COS. 



ARL1 resistance to other antifungal compounds 
or stresses 

To test whether overexpression of ARL1 confers resist- 
ance to other antifungal agents, cationic compounds, or 
osmotic stresses, wild type (vector control) and ARL1 
overexpression strain cells were grown in the presence 
of these perturbations. Overexpression of ARL1 did not 
confer resistance to NaCl, sorbitol or LiCl compared 
with the vector control (data not shown). Overexpres- 
sion of ARL1 does confer resistance to COS-5.44 at high 
concentrations (112.5 ug/ml) where the vector control is 
unable to grow (Figure 7A). Interestingly, the ARL1 
overexpression strain is as sensitive to Amphotericin B 
(12 ug/ml, Figure 7B) and Terbinafine (16 (not shown), 
8 and 4 ug/ml, Figure 7E and F) as the wild type cells. 
The ARL1 overexpression strain appears to be slightly 
more sensitive than the vector control to Fluconazole 
(32, 28, (not shown) 24 and 20 ug/ml, Figures 7C and 
D). These results suggest that the resistance to COS 
resulting from ARL1 overexpression is specific and does 
not extend to other antifungal compounds that disrupt 
fungal cell membranes [40]. Overexpression of ARL1 
allows cells to withstand H2O2 perturbation slightly bet- 
ter (12%) than control cells suggesting that ARL1 over- 
expression may be providing some protection against 
oxidative stress (Figure 7G). 

Synergistic interaction of COS-5.44 with fluconazole 

Because COS-5.44 and Fluconazole have different modes 
of action, we examined the effect of treating cells with 
both compounds simultaneously to test if these two 
compounds might interact synergistically. We quantified 
the degree of interaction across different concentrations 
in a dose response matrix for each drug interacting 
with itself and with the other drug (Additional file 2: 
Table S8 - Sll). At concentrations where either drug 
only weakly inhibits cells growth (i.e. COS inhibits 
growth < 15% while Fluconazole inhibits < 9%) when used 
in combination, a dramatic decrease in cell growth is 
observed (Figure 8A,B). Using the area under the growth 
curve as a metric, the combination of COS-5.44 and 
Fluconazole act in a synergistic fashion since the 
combination of both drugs inhibits growth much more 
(i.e. 45 to 86% inhibition) than the sum of growth inhib- 
ition caused by each compound individually. 

Discussion 

We have applied multiple comprehensive chemoge- 
nomics assays in yeast to identify the gene targets of 
COS-5.44. Unlike other antifungal drugs tested with the 
same approach [24] where typically 3-10 genes are iden- 
tified as potential drug targets, we found a much larger 
gene set (39) of yeast deletion strains highly sensitive to 
COS-5.44. The most sensitive heterozygous deletion 
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strain in the pool will often be the drug target [24,26]. In 
our case, we identified 17 heterozygous deletion strains, 
suggesting that COS does not target a single specific 
protein as its antifungal mode of action. In a previous 
global fitness assay similar to the HIP-HOP assay, 101 
homozygous strains were sensitive after 9 h (~ 6-7 gen- 
erations) of exposure to 25 (ig/ml of 10 kDa chitosan. 
Very few heterozygous strains were identified as sensi- 
tive to chitosan [23]. The differences seen in our results 
compared to those of Zakrzewska et al. [23] are likely 
due to a number of different experimental parameters. 
First, there were differences in the form and concentra- 
tion of COS used, (i.e. a lower concentration and larger 
molecular weight of COS was used in that study [23]). 



As we described earlier, the physiochemical properties of 
COS such as degree of deacetylation and molecular 
weight are known to affect COS's biological activity 
[16,41,42]. Second, the treatment period used in the 
Zakrzewska et al. study was shorter than in the one used 
in our study (9 vs. 40 hrs equivalent to ~ 6-7 and 20 
generations respectively). The shorter exposure time 
could explain the low number of chitosan sensitive 
heterozygous deletion strains found by these authors 
since the subtle phenotypes of the heterozygous dele- 
tions strains often require longer growth periods to re- 
solve fitness differences between strains [43]. Finally, the 
lower threshold for the identification of sensitive strains 
(i.e. those having a log2 ratio of > 1.585) used by 
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Figure 9 Potential mechanisms of COS resistance resulting from the overexpression of specific genes identified in this study. In this 
model, ERG24 overexpression would reduce membrane permeability by altering the composition and fluidity of the plasma membrane. /W/J, 
BCK2, MSGS and RBA50 either alone or in combination, could detect COS binding to the membrane and subsequently induce transcriptional and 
additional changes that reduce membrane permeability. Overexpressing these strains that are resistant to COS results in a down regulation of 
transcripts related to energy production (e.g. ATP synthesis and mitochondrial activity; examples of specific down-regulated genes are listed in 
the blue boxes). A large number of genes involved in cell division (mitosis and meiosis), cell cycle and cytoskeleton are transcriptionally 
up-regulated in the COS resistant overexpression strains (examples listed in the orange boxes). Integrin refers to yeast proteins that have 
properties similar to a mammalian integrins. 
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Zakrzewska et al. 2007, [23] accounts for the higher 
number of homozygous strains they determined to be 
chitosan sensitive. 

Several of the genes sensitive to COS in our HIP-HOP 
assay and resistant in the multicopy suppression assay 
were enriched for proteins targeted to membranes. This 
is not surprising given that COS likely perturbs mem- 
brane integrity. COS might also be creating oxidative 
stress leading to the accumulation of damaged proteins 
that are degraded by the proteosome. Protein degrad- 
ation/proteosome functions were among the biological 
processes enriched in the HIP-HOP assays (Figure 2C). 

We selected 5 of the confirmed overexpressing strains 
that provide resistance to COS to uncover any changes 
in their gene expression profiles that might give insights 
into COS's mechanism of action. ARL1 (YBR164C) was 
found to be sensitive to COS as a deletion strain in the 
HOP assay and resistant in the MSP assay in this study 
and was found to be sensitive to chitosan as a homozy- 
gous deletion strain in a previous study [23]. We con- 
firmed that an ARL1 deletion strain was sensitive to 
COS-5.44 and overexpression of ARL1 conferred resist- 
ance to COS-5.44 (Figures 3B and 3D). We believe the 
identification of this gene does provide information into 
the molecular mechanism of COS-5.44. Aril is a G pro- 
tein and soluble GTPase and is a member of the Ras 
superfamily [28,44]. Aril is highly conserved in all 
eukaryotes with 65% homology to human ADP- 
ribosylation factor-like protein 1 (Aril). Aril has been 
shown to be associated with the tra«s-Golgi and is 
thought to be required for endosome-Golgi trafficking 
[45,46]. Arll-GTP recruits specific receptor proteins to 
the membrane surface by binding to their GRIP 
domains, although in yeast, no specific Aril binding re- 
ceptor proteins have yet been identified [44,45]. Muta- 
tions in yeast ARL1 are not lethal [28] but ARL1 
mutants do show mild defects in localizing proteins to 
vacuoles as well as defects in potassium uptake [44,46-48] . 
aril mutants have also been shown to be more sensitive 
to antifungal agents such as Hygromycin B [46] . Aril can 
be myristolayted and we hypothesize that through myris- 
toylation, the soluble form of the GTPase could be bound 
to the membrane bilayer. In this way ARL1 could act as a 
sensor and modulate membrane trafficking at the onset of 
COS 5.44 induced membrane permeabilization (Figure 9). 
Aril (Arllp, PDB id lmoz) has been predicted to weakly 
associate with membranes (deltaG of -4 kcal/mol without 
ligands for the dimer) with the N-terminal residues (N- 
myristoyl glycine) of each dimer binding to the membrane 
[27,49,50]. GTPases are known to be involved in signal 
transduction pathways in filamentous fungi and yeast. 
These proteins are also involved in conidiation, a process 
that is enhanced by chitosan treatment in filamentous 
fungi [8]. When the ARL1 gene is overexpressed and yeast 



cells are challenged with COS, a large set of genes related 
to cytoskeleton organisation (e.g. microtubule dynamics) 
and stress sensing are up-regulated (Figure 5). These 
genes are also involved in cell division (mitosis and mei- 
osis) and the cell cycle. This could explain the ability of 
chitosan to enhance sporulation [28,44]. 

ERG24, which encodes C-14 sterol reductase, is a gene 
involved in ergosterol biosynthesis. Because ergosterol is 
a unique lipid in fungal membranes it is a classical target 
of many antifungal compounds. In a previous study, 
while ergosterol content was not associated with resist- 
ance to chitosan, in some filamentous fungi saturation of 
free fatty acids was associated with resistance [51]. Free 
fatty acid desaturases lower membrane fluidity and desa- 
turase mutants have been found to confer resistance to 
chitosan in N. crassa [51]. Membrane patches rich in er- 
gosterol and other lipids alter membrane fluidity and are 
important for correct membrane function [52]. It also 
has been shown that an increase in sterol biosynthetic 
capacity results in an increase in the availability of fatty 
acids [53]. We can envision that overexpression of 
ERG24, a key gene in the ergosterol biosynthetic path- 
way, could increase the size/presence of ergosterol rich 
membrane rafts and fatty acid content making the mem- 
brane less fluid and therefore more resistant to COS. 

For both yeast and filamentous fungi, the plasma 
membrane is the primary target for chitosan. In yeast, 
chitosan causes cell leakage and stops cell growth [22]. 
In filamentous fungi, chitosan permeabilizes the mem- 
brane in an energy dependent manner and kills cells of 
sensitive species such as N. crassa conidia [54]. A recent 
study by Palma-Guerrero et al. 2010, [51] found that 
membrane fluidity is a key factor in the antifungal action 
of chitosan with fluid fungal membranes binding more 
chitosan than less fluid ones. Chitosan crosses the fungal 
cell wall and binds to membrane lipids with a preference 
for negatively charged ones [51]. Chitosan permeabilizes 
the membrane when the concentration reaches a critical 
level [55,56]. Our study identifies pathways involved in 
membrane structure and signalling that are affected by 
COS. The Sytox uptake experiments confirm that mem- 
brane permeability is increased when cells are exposed 
to COS and that overexpression of ARL1 suppresses the 
increase in membrane permeability caused by COS 
(Figure 6). 

While our work and previous studies have shown that 
chitosan induces stress responses in yeast and other 
fungi, our study shows that pre-treatment with any of 
the well characterized environmental stresses (thermal, 
salt, osmotic and oxidative) does not provide resistance 
to COS-5.44 in yeast (Additional file 1: Figure Sll). 
Therefore, the changes induced by COS would appear to 
be different from previously described environmental 
stress response pathways. 
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Assuming that COS's main mode of anti-fungal action 
is to increase the permeabilization of membranes, we 
speculate that several of the genes we have identified in 
this study are involved in a cellular response that 
reduces permeabilization of the membrane and/or miti- 
gates the damage induced by the permeabilization. For 
example, ARL1, BCK2 and MSG5 could be part of a 
sensing and signalling mechanism(s) that detects COS 
binding to the membrane that induces transcriptional 
and other changes to reduce membrane permeability. 
ERG24 could be reducing membrane permeability by 
affecting the composition and fluidity of the membrane 
(Figure 9). Ultimately, if COS concentrations are too 
high, the cellular response would be insufficient to pre- 
vent membrane permeabilization and cell death would 
ensue - possibly as a result of the release/leakage of oxy- 
gen radicals from the mitochondria. While the exact 
mechanism of cell death by COS is still unknown, it 
appears to be a respiration/ ATP-dependent process in 
filamentous fungi [54]. Overexpressing strains resistant 
to COS display down-regulated genes related to energy 
dependent processes (ATP synthesis, mitochondrial ac- 
tivity, etc.) (Figure 5; Additional file 1: Figures S5-S9; 
Additional file 2: Table S6). In this sense, the overexpres- 
sing strains would mimic the effect of treatments such 
as sodium azide and low temperature that inhibit oxida- 
tive respiration and chitosan damage to fungal cells [54]. 
Interestingly, the ARL1 overexpression strain also shows 
increased resistance to certain forms of oxidative stress 
such as H 2 0 2 (Figure 7G). The transcriptome analysis of 
several of the overexpressing strains suggests that resist- 
ance to COS involves generation of less ROS by redu- 
cing oxidative respiration as well as reducing cell 
responses to ROS. 

A comparison of the gene targets found in the HIP- 
HOP screens in this study to genes found in similar 
screens for other fungicides that act at the cell mem- 
brane (e.g. Amphotericin B, Fluconazole and Terbina- 
fine) suggests that COS may have a different mechanism 
of action than these other drugs [57]. In support of this 
hypothesis, we have found that the ARL1 overexpression 
strain does not provide resistance against Amphotericin 
B, Fluconazole or Terbinafine (Figure 7B-F). This sug- 
gests that COS might be an effective antifungal agent in 
strains that are resistant to the other perturbing cell 
membrane compounds. 

Fluconazole, Terbinafine and Amphotericin B antifun- 
gal activity involves the targeting of ergosterol, the prin- 
cipal cell membrane sterol of fungi, by different 
mechanisms. Fluconazole (an triazole) inhibits 14a- 
demethylase (lanosterol demethylase, Ergll), a fungal 
cytochrome P450 dependent enzyme, and depletes cell 
membrane ergosterol resulting in reduced membrane 
fluidity, and the accumulation of 14a-demethylated 



sterols that are toxic to the cell. This leads to growth ar- 
rest and eventual fungal cell death [58,59]. Terbinafine 
inhibits squalene monoxygenase (Ergl), an enzyme re- 
sponsible for conversion of squalene epoxide, a precur- 
sor to lanosterol in the ergosterol synthesis pathway 
[59]. Amphotericin B targets the fungal cell membrane 
by directly binding to ergosterol, forming complexes that 
intercalate the cell membrane, resulting in formation of 
pores and leakage of intracellular contents [60,61]. COS 
appears to perturb the cell membrane possibly through 
the ergosterol pathway, by targeting ERG24, a gene that 
encodes for C-14 sterol reductase [62]. 

Because COS seems to have a distinct mode of ac- 
tion compared with other antifungals, we were inter- 
ested in seeing its inhibitory effect when used in 
conjunction with another antifungal drug. We tested 
COS and Fluconazole in combination (as well as on 
their own), and saw a strong synergism between them. 
That is, at concentrations where either drug alone had 
little inhibitory effect, a strong inhibitory effect that 
was much greater than additive was observed when 
they were combined (45 to 86% growth inhibition; 
Figure 8A,B). One possible explanation for this synergy 
might be that the increased membrane permeability 
induced by COS allows more Fluconazole to penetrate 
the cell membrane to inhibit 14a-demethylase that in 
turn further reduces membrane fluidity and increases 
the production of toxic sterols. 

Conclusions 

The results of this study have provided insights into the 
molecular mechanisms of the antifungal modes of action 
of COS. Some of the pathways COS appears to affect in- 
clude sensing, signalling, and the composition of the cell 
membrane (Figure 9). Our results suggest that COS does 
not have a single specific gene target. We provide evi- 
dence that COS increases cell membrane permeability 
and that the overexpression of ARL1 can reduce the per- 
meability induced by COS. Resistance to COS does not 
confer resistance to other antifungal compounds that 
disrupt the fungal cell membrane such as azoles, poly- 
enes and Terbinafine. COS can synergize with other 
antifungal agents such as Fluconazole suggesting that 
COS could be considered as alternative antifungal treat- 
ment to fungal pathogens that are resistant to other anti- 
fungal compounds or be used in combination with other 
antifungal agents to enhance their activity. 

Methods 

Chitosan oligosaccharide 

Chitosan oligosaccharide (COS-5.44) (5.44 KDa, PDI (poly 
dispersity index) = 1.14 and 97% degree of deacetylation) 
was kindly provided by Dr V. Tikhonov (Laboratory of 
Physiologically Active Biopolymers, A. N. Nesmeyanov 
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Institute of Organoelement Compounds, Russian Acad- 
emy of Sciences, Moscow, Russia). COS-5.44 was pre- 
pared as described by Tikhonov et al. 2006, [20]. For the 
experiments performed in this study, a COS-5.44 stock 
solution (25 mg/ml) in 50% DMSO was prepared, ali- 
quoted and stored at -20°C. 

Chitosan 

Chitosan (T8s) with a molecular weight of 70 kDa and 
80% degree of deacetylation was obtained from Marine 
BioProducts GmbH (Bremerhaven, Germany). Chitosan 
was prepared as described by Palma-Guerrero et al. 2008, 
[55]. Each batch of dialyzed and autoclaved chitosan was 
stored at 4°C for a maximum of 2 weeks. 

Yeast sensitivity to COS-5.44 and chitosan 

A pre-screen of COS-5.44 and chitosan (T8s) in wild 
type yeast (BY4743) was performed as described in 
Pierce et al. 2007, [43] to determine the antifungal activ- 
ities of these two compounds that have different mo- 
lecular weights and degrees of deacetylation. For the 
pre-screen of COS-5.44, concentrations from 62 ug/ml 
to 125 ug/ml of COS-5.44 were tested (Figure 2B). 98 ul 
of cell mix at an initial OD 600 of 0.0625 in 0.5X YPD pH 
5 (5 g bacto yeast extract, 10 bacto peptone and 10 g 
dextrose in 1 1 distilled water [63]) was added to each 
well of a 96 well plate (Nunc) then 2 ul of the corre- 
sponding COS-5.44 stock were added, including carrier 
controls (50% DMSO and 0.25 mol/1 HC1 neutralized 
pH5.6). The chitosan (T8s) screen followed the same 
protocol, with the slight difference in volume of cell mix 
(97 ul) added and chitosan (T8s) stock (3 ul). Con- 
centrations from 125 |ig/ml to 162.5 ug/ml of chitosan 
(T8s) were tested (Additional file 1: Figure SI) Cells 
were grown at 30°C in a shaking-spectrophotometer 
(Tecan), with readings taken every 15 min for 24 hrs. 
Tecan ODs were converted to conventional 1 mm path 
length cuvette ODs using a calibration function provided 
by Ericson et al. 2010, [24]. The COS-5.44 pre-screen 
was performed to identify an IC 70 _9o and ICi 0 _20 of 
COS-5.44 in the BY4743 to be used in the multicopy 
suppression and deletion profiling (HIP-HOP) respect- 
ively (Figure 2B). 

Hapioinsufficiency and homozygous profiling 
(HIP-HOP) assay 

Hapioinsufficiency and Homozygous profile assays, gen- 
omic DNA purification, PCR amplification of barcodes, 
array hybridization, analysis and confirmation were per- 
formed as described in Ericson et al. 2010, [24]. The 
homozygous collection was grown for 5 generations (ap- 
proximately 10 hrs) and the heterozygous collection for 
20 generations (approximately 40 hrs), as these periods 



of time have been found to be optimal to allow each col- 
lection to resolve growth differences [43]. 

Experiments were performed in duplicate. Yeast dele- 
tion strains with log2 ratio of 3.5 or higher (Highly sen- 
sitive COS-5.44 deletion strains) were selected for 
individual confirmation. A log2 ratio of 3.5 represents an 
approximately 11- fold less abundance of the deletion 
strain after COS treatment. 

Multicopy suppression profiling and screen conditions 

MSP screen, plasmid isolation, Affymetrix TAG4 ar- 
ray hybridization, microarray results analysis and con- 
firmation followed the protocols described by Hoon 
et al. 2008, [25] and Ericson et al. 2010, [24]. 

Transformation of candidate resistant genes for 
overexpression in yeast 

Sixty-eight genes were identified with MSP as possible 
candidates to confer resistance to COS-5.44 (Additional 
file 1: Figure S2). To confirm which of the 68 genes con- 
ferred resistance, each gene was individually overex- 
pressed in yeast. For this purpose E. coli strains carrying 
the genes of interest were obtained from the MoBY-ORF 
2 u collection, kindly provided by Sarah Barker from the 
Boone lab [30]. 57 of the 68 candidate genes were avail- 
able in the MoBY-ORF 2 u collection and thus 11 of the 
genes of interest were not present in this collection and 
could not be tested. 

Individual colonies of each E. coli strain were grown in 
selective media in 2YT (1% Yeast extract, 1.6% Bacto 
tryptone, and 0.5% Sodium) + 0.4% Glucose + Carbenicil- 
lin (200 mg/ml, SIGMA) + Kanamycin (50 mg/ml, 
SIGMA) in 96 square well blocks (Greiner, Bio-one) at 
37°C and 200-250 rpm. Cells were harvested (5 min, 
4000 rpm) to proceed with DNA extraction. Macherey- 
Nagel Multi-96 purification kit (MN Cat # 740 625.4) 
was used to miniprep the plasmids, following the manu- 
facturer's protocol. Elution was done in 150 ul to 
collect ~ 100 ul of purified DNA. 

Each plasmid extract was double digested with EcoRI/ 
Xhol (New England BioLabs Inc.). Digested samples 
were run on a 0.8% agarose gel to confirm vector and 
ORF fragments size. 

The confirmed plasmids were transformed into yeast 
strain BY4743 (his3 l/his3 1 leu2 0/leu2 0 ura3 0/ura3 0 
metl5 0/+ lys2 0/+), using a 96 well yeast transformation 
protocol as described by Gietz et al. 2007, [64]. Five 
single colonies of each strain were grown overnight 
in selective liquid media then archived in 15% glycerol 
at -80°C. 

Once overexpressing strains were obtained, confirm- 
ation of resistance to COS-5.44 was performed 
(Figure 1). Overexpressing strains were grown at differ- 
ent concentrations of COS-5.44 previously found to be 
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inhibitory to the wild type strain in MSP. Concentrations 
tested were 101.25, 106.28, 112.5, 118.75, 125.0 and 
250.0 ug/ml. To ensure that the resistance was conferred 
by the overexpressed gene, plasmid DNA was extracted 
from two colonies of each overexpressing strains, and 
transformed into E. coli. Amplified plasmids were 
extracted using Genejet kit (Fermentas) following manu- 
facturer protocol. Purified plasmids were double digested 
as described. Confirmed plasmids were retransformed 
into BY4743 and three single colonies were selected 
and confirmed for COS-5.44 resistance as previously 
described. 

Sample collection and RNA extraction of 
overexpressing strains 

A confirmed colony of each overexpressing strain and 
the wild type with the empty vector was grown over- 
night at 30°C and shaking at 200-250 rpm. An overnight 
culture was used to inoculate 3 flasks each with 400 ml 
of 0.5X YPD pH 5. Cells were grown until the culture 
reached an OD 600 of 0.8 and a 100 ml sample was col- 
lected (time zero). The remaining 300 ml were treated 
with COS-5.44 (112.5 ug/ml) and 100 ml samples were 
collected at 15, 30 and 60 min after treatment (Figure 1). 
Cells were collected by centrifugation (5 min, 4000 rpm) 
frozen with liquid nitrogen and stored at -80°C, for 
RNA extraction and labeling. Once all samples were col- 
lected a hot acid phenol-chloroform RNA extraction was 
performed [65]. 

Microarray analysis (see below) of the COS treated 
wild type cells found that the maximum number of dif- 
ferentially expressed genes were obtained with the 60 
min treatment with minimal loss of genes differentially 
expressed at earlier time points but not at 60 min. 
Therefore, for the overexpression strains, microarray 
analysis was only performed on the untreated and 60 
min COS treated samples. 

cDNA and microarray sample preparation 

cDNA synthesis and sample labeling for NimbleGen 
4X72k yeast microarrays (Roche NimbleGen, Inc.; De- 
sign ID A6186-00-01, TI4932 60mer expr X4) was per- 
formed as described in the manufacturer's protocol with 
minor modifications. cDNA synthesis was performed 
using 10 ug of total RNA. The sample labeling reactions 
were done using 1 ug of double stranded cDNA. Three 
cDNA biological replicates either with or without a 60 
min exposure to COS-5.44 for each of the 5 overexpres- 
sing strains as well as an untransformed wild type 
BY4743 cells (vector control; for a total of 36 samples) 
were hybridized to NimbleGen 4X72k microarrays 
(Roche NimbleGen, Inc.). Microarrays were scanned 
with a Genepix 4000B scanner (Molecular Devices Inc. 
Sunnyvale, CA). 



Microarray data analysis 

Raw images were obtained using GenePix Pro software 
(Version 5.0, Molecular Devices Inc. Sunnyvale, CA). Raw 
probe intensities (x, y and signal reports) were obtained 
using NimbleScan software (v2.4, Roche NimbleGen, Inc). 
Data preprocessing of 60mer oligonucleotide arrays was 
performed using the BioConductor package Oligo [66]. 
Genes with significant expression difference compared 
with the vector control (wild type, BY4743) were identified 
for each overexpressing strain using the BioConductor 
package Limma Linear models for microarray analysis, 
[67], taking P-value<0.05 and log2 fold change of > 1 or 
< -1 as significant unless otherwise noted. The data dis- 
cussed in this publication have been deposited in NCBI's 
Gene Expression Omnibus [68] and are accessible through 
GEO Series accession number GSE32888 (http://www. 
ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE32888). 

Hierarchical clustering 

Hierarchical clustering of the data was performed using 
the program Cluster 3.0 originally written by Eisen, [69], 
modified by de Hoon, [70]. Expression data of the 5 
overexpressing strains was clustered using city-block dis- 
tance measurement and average linkage clustering 
method. The resulting clusters were visualized using Java 
TreeView [71]. 

Functional annotation of chemogenomic assays and 
microarray results 

The genome-wide profile of the COS-5.44 sensitive dele- 
tion strains was examined using Gene Set Enrichment 
Analysis GSEA, [29] to identify enriched biological pro- 
cesses. GSEA has been considered advantageous for 
quantitative genome-wide profiles because it can take 
into account the entire profile, exploiting the weighting 
of the genes (e.g. degree of sensitivity to a drug (i.e. fit- 
ness defect), differential gene expression). Other func- 
tional analysis methods focus on a list of genes that 
satisfy the characteristic of interest. This would first re- 
quire setting a threshold to divide the profile into lists of 
genes that satisfy the characteristic of interest and those 
that do not. It is often unclear how to choose an appro- 
priate threshold, yet with GSEA, choosing such a thresh- 
old is not required. 

Deletion strains were mapped to genes using chromo- 
somal feature data downloaded from the Saccharomyces 
Genome Database (SGD) on April 16, 2011. Multidrug 
resistance genes MDR, [57] were filtered from the 
genome-wide profile of sensitive deletions strains to 
COS, to identify biological process specific to COS re- 
sponse. The list of all genes ranked by their fitness 
defects was analyzed by GSEA v2.07 (4864 genes, with- 
out MDR genes) [29]. Default parameters were used ex- 
cept that minimum and maximum gene set size were 



Jaime ef al. BMC Genomics 2012, 13:267 
http://www.biomedcentral.com/1471-2164/13/267 



Page 23 of 26 



restricted to 5 and 300, respectively. Biological process 
annotation was obtained from Gene Ontology website 
(http://berkeleybop.org/goose) on April 13, 2011. Add- 
itional protein complex annotations based on consensus 
across different studies were obtained from Benschop 
et al. 2010, [72]. 

Enrichment maps were generated with Enrichment 
Map Plugin vl.l [73] developed for Cytoscape [74] with 
default parameters. Nodes in the maps were clustered 
with the Markov clustering algorithm, using an overlap 
coefficient computed by the plugin as the similarity 
metric (coefficient < 0.5 were set to zero) and an inflation 
parameter with value of 2. For each cluster leading edge 
were computed as in Subramanian et al. 2005, [29] for 
each member of a node. The top 10 leading edge genes 
associated with the most nodes in the cluster are shown 
in the bar plot (if less than 10 genes, all are shown). 

A functional analysis (i.e. GSEA [29]) was performed as 
described above for the transcriptional profiles of the wild 
type (BY4743) exposed to COS-5.44 and the 60 min chito- 
san data set from Zakrzewska et al. 2005, [22] that was 
performed on X2180-1A, MAT* SUC2 mal gal2 CUP 
cells. The set of environmental stress response genes ESR, 
[75], were filtered from the transcriptome data sets to 
identify the biological processes enriched by COS or chit- 
osan treatment apart of the general stress response of 
yeast. Enrichment maps were generated with Enrichment 
Map Plugin vl.l [73] developed for Cytoscape [74] with 
the following parameters: FDR < 1%, p-value < 0.005 and 
using an overlap coefficient computed by the plugin as the 
similarity metric (coefficient < 0.5 were set to zero). En- 
richment maps were compared with the comparison func- 
tion of the Enrichment Map Plugin vl.l [73]. 

A similar functional analysis was performed on the 
whole transcriptome dataset of each overexpressing strain. 
As was described above, the set of environmental stress 
response genes were filtered from the transcriptome data 
set. An additional enrichment analysis of significantly up- 
or down-regulated transcripts (P-value > 0.05 and log2 fold 
changes >1 or<-l) was also performed using the single 
enrichment analysis tool of Babelomics [35,36]. All default 
parameters were used. 

Quantitative RT-PCR 

Quantitative reverse transcription - PCR (qRT-PCR) 
verification of the microarray results for selected genes, 
was performed using KAPA SYBR FAST (KAPA Biosys- 
tems) in a CFX384 real-time PCR detection system 
(Bio-Rad), with two biological replicates from wild type 
(vector control BY4743) or Aril overexpressing strain 
that were either untreated or exposed to COS-5.44 for 
60 min. Each RT qPCR reaction had 200 pg of cDNA, 5 
uL of KAPA SYBR FAST master mix (2x), 300 uM of 
each primer and water to have a final volume of 10 uL 



per reaction. Cycling conditions used were as follows: 
95°C for 3 min, followed by 40 cycles of 95°C for 3 s and 
60°C for 30 s then 95°C for 10 s. All the reactions were 
performed in triplicate, using primers for the following 
genes: ACT1, ARL1, CMD1, COX5B, ERV2S, FIG2, 
HSP30, MUC1, RCR1, RPL32, SPT15, UBI4 and VBA5. 
Primers were designed using Primer-BLAST Developed 
at NCBI uses Primer 3, [76]. See Additional file 2: Table 
S6 for the sequences of the primers used. The average of 
three technical replicates was normalized to four in- 
ternal control transcripts, ACT1, ERV2S, RPL32, SPT1S, 
there were not significant differences among internal 
controls therefore all were used for normalization. Data 
was analyzed using CFX Manager (Bio-Rad laboratories, 
Inc). 

Cell permeability assay 

Plasma membrane permeability was measured by 
SYTOX Green uptake as described by Thevissen et al. 
1999, [37] with modifications. S. cerevisiae cells of the 
wild type and overexpressing strain Aril were grown in 
the presence of COS-5.44 (91.1, 101.25, 106.28 and 
112.5 mg/ml) and vehicle control (1% DMSO). Cultures 
of the overexpressing strain Aril and wild type (B4743) 
with empty vector were set up at initial OD of 0.0625 
and grown overnight (17 hrs). After COS-5.44 treatment 
the cells were pelleted by centrifugation (5 min, 
4000 rpm) and washed three times in 0.1 M Tris-HCl, 
pH 7.0. Two hundred microlitre aliquots of the yeast cell 
suspension was incubated with 0.2 uM SYTOX Green in 
96-well microplates for 30 min at 30°C with periodic agi- 
tation in the dark. Three replicates of each strain were 
performed. Fluorescence was measured using a Guava 
easyCyte flow cytometer (Millipore). The green fluores- 
cence filter set 488 nm (excitation) and 520 nm, (emis- 
sion) was used. Fluorescence data was analyzed with the 
software Flowjo (V9.3.1, Tree Star, Inc.). 

Stress response and COS-5.44 resistance 

To verify if the activation of the environmental stress re- 
sponse (ESR) provides resistance to COS-5.44 treatment, 
the wild type strain (BY4743) was exposed to sub-lethal 
doses of primary stresses previously reported by Berry 
and Gasch 2008, [38] to provide resistance to more se- 
vere stresses. The yeast wild type strain was exposed to 
sub-lethal doses of stress (mild primary stress) then cells 
were exposed to chitosan oligosaccharide treatment (se- 
vere secondary stress). The primary mild stresses tested 
were: 0.25 M and 0.5 M NaCl, 1 M sorbitol, 0.003% and 
0.006% H 2 0 2 for 60 min, 37°C for 15 min and 30°C con- 
trol. Yeast cells were collected by centrifugation (5 min, 
4000 rpm) then exposed to the secondary stress treat- 
ment COS-5.44 for two hours (91.1, 101.25, 0.1068 and 
112.5 ug/ml). After the secondary stress, cells were 



Jaime ef al. BMC Genomics 2012, 13:267 
http://www.biomedcentral.com/1471-2164/13/267 



Page 24 of 26 



collected by centrifugation, transferred to YPD and their 
growth rates measured. 

Resistance to other antifungal and cationic compounds 

To ask if the overexpression of ARL1 conferred resist- 
ance to other cationic compounds, the wild type 
and ARL1 overexpression strain were grown in the 
presence of other cationic compounds (as described in 
section 3.2). 

A dose response assay was performed in the wild type 
and overexpression strain ARL1, as previously described 
screen protocol, using the following compounds: NaCl 
(0.5 and 1 M); sorbitol (1 and 2 M); LiCl (0.01 M); 
hygromycin B (25, 12, and 6 ug/ml); Amphotericin B 
(25, 12 and 6 ug/ml); Fluconazole (32, 28, 24 and 20 ug/ 
ml); Terbinafine (16, 8 and 4 ug/ml); H 2 0 2 (0.012% 
0.006%); 1% SDS (0.05, 0.025 and 0.012%). 

Interaction between COS and fluconazole 

For COS-5.44, a concentration range of 35-105 ug/ml 
and for Fluconazole, a concentration range 4-24 ug/ml 
was used. Stocks for each drug treatment were prepared 
in 12.5% DMSO resulting in a final concentration of 1% 
DMSO. Wild type yeast (BY4743) was grown overnight 
and diluted at an initial OD 600 of 0.0625 in 0.54X YPD 
pH 5. 92 ul of cell mix was aliquoted to each well of a 
96 well plate (Nunc), and 4 ul of the corresponding 
COS-5.44 and/or Fluconazole stock were added in a 
7 x 12 dose response matrix. Two dose response matri- 
ces were performed: self-self (COS-5.44 - COS-5.44 and 
Fluconazole - Fluconazole) and two replicates of the 
COS-5.44 and Fluconazole matrix. Cells were grown at 
30°C in a shaking-spectrophotometer (Tecan), with 
readings taken every 15 min for 24 hrs. We used the 
area under the growth curve as a metric to measure 
growth inhibition and drug interaction. 
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